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Abstract 

The well-known Regge- Wheeler equation describes the axial perturbations of 
Schwarzschild metric in the linear approximation. From a mathematical point of 
view it presents a particular case of the confluent Heun equation and can be solved 
exactly, due to recent mathematical developments. We present the basic properties 
of its general solution. A novel analytical approach and numerical techniques for 
study the boundary problems which correspond to quasi-normal modes of black 
holes and other simple models of compact objects are developed. 

1 Introduction 

In the linear approximation the axial perturbations of Schwarzschild metric of Keplerian 
mass M are described by the well-known Regge- Wheeler equation pQ: 

d 2 ®^ + {-d 2 x + v Sjl ) $ S)Z = o. (l.i) 

Here x = r + r s In (r/r s — 1) is the Regge- Wheeler "tortoise" coordinate, r s = 2M(= 1) 
is the Schwarzschild radius, and $> s ,i(t, r) is the radial function for perturbations of spin 
s and angular momentum I > s. The corresponding effective potential V s j(r) reads: 

V s j(r) — (l — i) (^—^ + ^s~j ■ The area radius r in it can be expressed explicitly as a 

function of the variable x using Lambert W-function [2\: r/r s = LambertW (e x ^ Ta ~ 1 ) + 1. 
The standard ansatz $ s ,i(t, r) = R £jS j{r)e %et brings us to the stationary problem: 

d?R £>S)l + (e 2 - V S)l ) R e>s>l = 0. (1.2) 

A large amount of references on this subject can be found in the review articles jSJ. 
There we meet often the statement that the exact solution of the problem is not known. 
Therefore, quite sophisticated approximate analytical and numerical techniques for solu- 
tion of physical problems related to Eqs. (jl.ljl . ()1.2j) were developed in the past. They 
discovered rich and important mathematical and physical properties of the corresponding 
problems. Some of them need further justification. 

A hint on the class of new mathematical functions, which solve the stationary Regge- 
Wheeler problem (jl.2j) is the use of the very time-time component of Schwarzschild metric 
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gu = 1 — 1/r = g as an independent variable, instead of the area-radius r. This transforms 
the Eq. JT2J) into: 



with the potential W a i(g) = l(l+l)g(l— g) 2 + (l — s 2 )g(l — g) 3 which describes anharmonic 
oscillators of different type, depending on the values of spin: s — - for scalar waves, 
s — 1 - for electromangnetic waves, s — 2 - for gravitational waves (see Fig. [!}. 



Figure 1: Typical forms of the potentials W a> i(g) for different s,l. The physical domain 
is g £ (0, 1). The point g = represents the Schwarzschild event horizon, the point g = 1 
- the physical infinity. 

This simple observation sheds additional light on the physical content of Eq. (|1.2|) and 
yields a better understanding of the qualitative differences and similarities between the 
perturbations of different spin. In particular, one can see in a more clear way the behavior 
of the potential, the RW equation and its solutions at infinite space points when these 
are placed at finite distances, using the variable g. 

Here we utilize the form (jl.3|) of the problem at hand only as a hint on a new direction 
of mathematical investigation. It is well known that the above polynomial form of the po- 
tential W 8 ,i{g) leads to a solution of geodesic equations in Schwarzschild metric in elliptic 
functions. According to the early article jl], in this case the corresponding Schrodinger 
like spectral problem can be solved by the method of continuous fractions. This method 
was developed for the specific RW spectral problem by Leaver [3| and at present is one of 
the basic methods for finding the QNM spectra. 

The modern mathematical development shows that many physical problems, including 
anharmonic oscillator ones, are solved exactly by Heun's functions [H]. Thus the form of 
the Eq. (|1.3|) suggests to describe the exact solution of Regge- Wheeler equation using 
Heun's functions. We shall show that the Eqs. (jl.2j) . (jl.3J) can be transformed to a 
particular case of confluent Heun equation. 

The confluent Heun equation of general form is well known and studied in some detail 
in recent mathematical literature j^j. At present its solutions are more or less well- 
described special functions, already included in modern computer packages. In particular, 
the Heun's functions are included in the last version of the widespread computer package 
Maple 10, which has been used in all original numerical calculations in our article. These 



d 2 x R £>s>l + (e 2 - W a>l ) R £ , s ,i = 



(1.3) 
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functions present a nontrivial generalization of familiar hypergeometric functions and have 
much richer and complicated properties, since the Heun equation has one more singular 
point, than the hypergeometric equation. 

In the present article, we give the exact solution of the Regge- Wheeler equation in 
terms of Heun's functions and apply them for study of different boundary problems in 
the relativistic theory of gravity. 



2 Solutions of Regge- Wheeler Equation 

2.1 Local Solutions of Confluent Heun Equation 

It is convenient to use the standard variable r for reaching the standard canonical form 
of the Heun equation of our problem. The explicit form: 

d 2 R £ ,s,i | 1 dR e>8il | e 2 r A -l(l + l)r 2 + (l(l + l) + s 2 -l)r + l-s 2 R ; = ()(21) 
dr 2 r(r — 1) dr r 2 (r — I) 2 £ ' s ' 

of Eq. (jl.2j) shows that it has three singular points: the origin r = 0, the horizon r = 1 
and the infinite point r = oo [Sj. The first two are regular singular points and may be 
treated on equal terms. The last one is an irregular singular point, obtained as a result 
of confluence of two regular singular points in the general Heun equation. Hence, we 
are to consider a specific representative of the class of confluent Heun equations. This 
circumstance is the source of complicated analytical properties of solutions of the Regge- 
Wheeler equation, as we shall see below. 

For transformation of Eq. (J2.1|) to a canonical form of the Heun equation j3] : 

d 2 H ( f3 + l 1 + l\dH 
+ [a+- + - — - — + 



r(r 



dr 2 \ r r — 1 J dr 

^ T) (^(5 + la((3 + 1 + 2)y + V +^ + ^( 1 -a)(P + l)^H = (2.2) 



one may use the ansatz R EtS j(r) = r s+1 (r — l) l£ e l£V H(r). It yields the following specific 
parameters: a = 2ie, (3 = 2s, 7 = 2ie, 5 = 2e 2 , 77 = s 2 — 1(1 + 1). 

There exist several substitutions of that kind. All they transform Eq. (j2.1J) to the same 
form (J2.2)) . but with different coefficients, for example, with all possible combinations of 
the signs of the coefficients a, f3 and 7. The corresponding solutions of equations (|2.2|) . 
with different coefficients, are related by simple transformations, described in literature 
jS]. It is useful to know all sets of such transformations and corresponding parameters. 
In some cases this permits us to construct two independent solutions of stationary Regge- 
Wheeler equation (jl.2j) . using the basic standard canonical solution of confluent Heun 
equation IfOjl 1 : 

HeunC (a, f3,j,8,r),r) , (2.3) 

which is defined by the convergent Taylor series expansion around the origin r = and 
present article we apply the notations for Heun's functions, which are used in Maple 10. 



normalization: 



HeunC (a, (3, 7, 5, t], 0) = 1, (2.4a) 
■^HeunC (a, /?, 7, 5, rj, 0) = -j-^ '-. (2.4b) 

The Taylor series expansion of the standard canonical solution ()2.Hj) with respect to 
the independent variable r is obtained using the known three-terms recurrence relation 
|.5J and initial conditions (J2.4|h 

It is not excluded that some specific sets of parameters a,/3, 7, 5, 77 may not yield 
well defined canonical solutions ()2.3|) of the confluent Heun equation ()2.2|) . This indeed 
happens in our problem, since in it s G Z and s > 0. As a result, for example, solutions 
HeunC (a, (3, 7, 5, n, r) with (3 = —2s are not well defined, because the denominator of the 
coefficients in Taylor series expansion includes products of form /3(/3+l)...(/3+N), N G Z + . 

One can go outside the domain of convergence of the Taylor series expansion - the 
open unit disk in C r , centered at the origin, using an analytical continuation in this 
complex plain. The relations between Heun functions of different values of parameters 
and arguments, described in the mathematical literature [^j, are useful for this purpose, 
too. 

Using these techniques we obtain the following local solutions of the confluent Heun 
equation ()2.2|) : 

1. Around the singular point r = there exist two independent local solutions H^^(r) 
with asymptotics: 

# ( 2>)~-L n®-{r)~l. (2.5) 

The second one is the standard canonical solution: 

H { °l~(r) = HeunC (2ie, 2s, 2ie, 2e 2 , s 2 - 1(1 + 1), r) , (2.6) 

which has convergent Taylor series expansion around the origin r = 0. This solution is 
the only one, which is finite at the point r = 0. It is normalized according to the relations: 



W°> s l (0) = HeunC (2ie, 2s, 2ie, 2e 2 , s 2 - 1(1 + 1), 0) = 1, (2.7a) 
Atf ( £ °)-(0) = ^HeunC (2ie, 2s, 2%e, 2e 2 , s 2 - 1(1 + 1), 0) = + + (2.7b) 

In the general case, the solutions of this type diverge at the horizon r = 1. There 
exist a countable number of such solutions which are finite simultaneously at the origin 
and at the horizon. Considering them, one can derive a discrete spectrum of values E n ^ s j 
for the boundary problem on the interval r G [0, 1]. To the best of our knowledge, this 
boundary problem has not attracted attention in the existing physical applications and 
its spectrum is not known at present, although it must be interesting for the study of 
the interior of Schwarzschild black holes 2 . Indeed, these solutions are obviously square 
integrable functions on the interval r G [0, 1]. If they will be proven to form a complete 

2 Note that in this domain the variable r plays the role of a time variable and the dependence of the 
former stationary solutions on the new space variable t is described by the simple exponential factor e let . 



set in the corresponding functional space, these solutions will present a natural basis 
of normal modes for a well defined expansion of any small perturbation of the metric 
in the interior of the black hole. Therefore these solutions seem to deserve a separate 
consideration. 

The first independent local solution in Eqs. (|2.5jl - ^(r) cannot be described by 
standard solution (|2.3j) of the confluent Heun equation for any set of parameters. Since 
b = 2s is positive integer, the second solution includes the In r terms [3] . It is analogous to 
the corresponding case of solutions of confluent hypergeometric equation, but we were not 
able to find more information about this solution in the available mathematical literature. 
Taking into account the well-known example of hydrogen atom, described by confluent 
hypergeometric functions, one can think that this solution may not play an essential role 
in most of the physical applications. Nevertheless, for completeness we shall add formally 
this solution to our list, using for it the notation H^^r). 

2. We have two independent local solutions of the Heun confluent equation ()2.2|) : 

#W+(r) = HeunC (-2ie, +2ie, 2s, -2e 2 , 2e 2 + s 2 - 1(1 + 1), 1 - r) , (2.8a) 
H { V~(r) = (r- l)- 2i£ HeunC (-2ie, -2ie, 2s, -2e 2 , 2e 2 + s 2 - 1(1 + 1), 1 - r) . (2.8b) 

The functions if*^(r) and (r — l) 2m H ^ s ; (r) have convergent Taylor series expansion 
around the second singular point - the horizon r = 1. They are normalized in the 
following way: 



{r-lf K H%(r) 
f(( r -l)^H%{r] 



2e 2 (l =f 1) ± 2ies+ie(l ± l) + s(s + l)-l(l + l) 
±2ie + 1 ' 



(2.9a) 



(2.9b) 



The first solution (J2.8a|) is the only one which is regular at the horizon r = 1. In 
general case, the two solutions (|2.8|) diverge at both the origin r = and the infinity 
r = oo. Using them one can formulate the corresponding two-point boundary problems 
and find their discrete spectra. We shall consider some of them in the next sections. 

3. Using known three-term recurrence relations for the coefficients a^, one is able 
to construct two asymptotics-series solutions around the infinite irregular singular point 
r = oo jHj: 



~ r- s - l e- 2l < r+ ^ J2^7> a ° = 1 ( 2 - 10a ) 
H ( S,7(r)~r-^J2% fl o=l;- (2.10b) 



u>0 



They are valid in the complex domain 



- | - /(/ + 1) + s 2 < arg(r) + arg( £ ) < 3 1 + 1(1 + 1) - s 2 , (2.11) 

see the articles by A. Decarreau et al. in [S], where one can find the restrictions of this kind 
for the general case of the confluent Heun equation. In our particular problem Eq. (J2.11|) 



leads to inessential restrictions on the variables in the complex plain C, since the length of 
the admissible interval of arguments arg(r) and arg(e) is 27r+2/(/ + l)— 2s 2 > 2(7r+s) > 2n 
for I > s > 0. 

Unfortunately, the transition coefficients between the three sets of local solutions (|2.6|h 
(12. 8j) and (|2.1()j) are not known at present. The only available information is negative: 
these coefficients cannot be expressed even in terms of hypergeometric functions [Sj. The 
same transition coefficients relate the corresponding local solutions of the Regge- Wheeler 
equation, see below. 



2.2 Local Solutions of Regge- Wheeler Equation 

In accord with the above local solutions of the Heun equation (12. 2j) . we obtain the following 
three sets of independent local solutions of the Regge- Wheeler equation (jl.lj) : 

^ > ,t l (t,r) = e"'R%(r), (2.12a) 

/<>> 



i2 (0)± O) = r s + 1 e ier+i£ln(1_r) < £ ' s ' 1 (2 12b) 

£ ' sA ' [HeunC(2te,2s,2te,2e 2 ,s 2 -1(1 + l),r); V ' ; 



*%(t,r) = e i * t R%(r), (2.13a) 
#W±(r)=r s+ y er±ielnM (2.13b) 



<J(^ = e^>), (2.14a) 
R { ?}f{r) ~ — ^ 4 = 1- (2.14b) 

They are correspondingly defined in domains around the three singular points of the 
problem: the origin r = 0, horizon r = 1 and infinity r = oo, as pointed in the upper 
indices of our notation. Any other radial function R £ yS ,i{ r )-, being solution of the sta- 
tionary Regge- Wheeler equation (jl.2)l . can be represented locally in the form of a linear 
combination: 

Re, s ,i(r) = C x+ Rf+(r) + C x ^-(r), (2.15) 

where Cx± are proper constants and X — 0, l,oo indicates the corresponding singular 
point. 

For fixed values of the parameters e, s, I each of the three sets ()2.12|) . (|2.13J) and ()2.14|) 
includes two independent local solutions of the stationary Regge- Wheeler equation (jl.2|) . 
Hence, in the corresponding common domains of validity of these solutions one obtains: 

R [ *S(r) = (e, *, + r£(e, a, /)<],7(r), (2.16) 

expanding one set of local solutions onto another. The main obstacle for analytical treat- 
ment of different interesting problems related to the Regge- Wheeler equation is that the 
transition coefficients fy^e, s, I) for X, Y = 0, 1, oo; I/F are unknown explicitly. 



2.3 Some Basic Properties of the Local Solutions of Regge- 
Wheeler Equation 

2.3.1 In-Out Properties of the Local Solutions 

Taking into account the relation (|2.4a|) . one easily obtains directions of spreading of 
the corresponding waves described by local solutions (|2.12a|) . ()2.13aj) and ()2.14aj) in the 
vicinity of the corresponding singular points X when the time t increases 3 . These "in" 
and "out" properties of local solutions are illustrated by the graphical schemes, shown in 
Tabled 



<5(*,r): 


«- 




-> 




1 <- 




1 -> 




— > oo 




<— oo 



Table 1: In-Out properties of local solutions 



Thus, we see that according to our conventions, the local solutions $ e j z describe 

waves going in the singular point X and the local solutions $ es i describe waves going 
out of singular point X. 

2.3.2 Time and Space Limits of Local Solutions for Complex Values of e 

In general, we have to consider complex valued e = Sr + iei G C £ . Since in the Regge- 
Wheeler equation only the squared quantity e emerges, without loss of generality one can 
accept conventionally that its real part is nonnegative: Er > 4 . 

The local solutions (|2.12af . ()2.13a|) and (j2.14aj) have the following obvious time-limits: 

( Vr, if e I > and |i? ( f s ^(r)| ^ oo; 

^ 1$ ( 2,f (t,r)| = |<ij(r)|, if e/ = 0; (2.17) 

[oo Vr, if ej < and \R { fJf(r)\ ^ 0. 

Hence, from a physical point of view we have to work only with ej > PJ, because 
the solutions with £j < are certainly unstable and, therefore, nonphysical. Then, for 
finite values of t in general we obtain the following space limits of the local solutions: 

lim|$ ( Jj-(t,r)| = oo, (2.18a) 
lim|$ ( Jj(t,r)| = 0. (2.18b) 

These need some justification for different singular points: 



3 The only still unstudied solution H ^ s + ;(r) has to be chosen in a proper way to insure the fulfilment 

of our "in-out" conventions. This requirement fixes the solution iJ^ s + ; (r) up to inessential normalization. 
4 In the case er = 0, we may have "pseudo" waves. The corresponding solutions do not oscillate. 



1. The limits (|2.18jl of local solutions around the origin r = are valid for s>0 and on 
any contour to this point in the complex plane C r . For s = the limit (j2.18b|) equals 1. 

2. Around the horizon r = 1 the limits (j2.18|) of local solutions are valid for all real 
values of s > and on any contour to this point in the complex plane C r . 

3. The behavior of limits (j2.18|) of local solutions in the complex plane C r around the 
infinite irregular singular point r = oo is more complicated and depends on the direction of 
the contour which approaches this point. This reflects the well-known Stocks phenomenon 
(see references in [HE]), and is related with the multiplier 

e ^fier-ielnr _ e i(a=F|e||r| cos(arg(r)+arg(e))) e ±|e||r| sin(arg(r)+arg(e))+|e| In |r| cos(arg(e)) 

in (|2.14b|) . Thus we see that 

j^Mi^l ^ e ±|e||r|sin(arg(r)+arg(e))+|£|ln|r|cos(arg(£))^ (2.19) 

and in the case |e| > we obtain: 

lim \R^Ur)\ = I °°' if arg(r) + arg(£) E (°^/^ od ^ (2 . 2 o a ) 
^" if arg(r) + arg(e) G (-7r,0)/mod 2tt. V ; 

if arg(r) + arg(e) G (0,7r)/mod 27r; 
if arg(r) + arg(e) G (— n, 0)/mod 27T. 



-CO 



+ 0O 




> \R { ^j(r)\ = c " 6 yj : ~*?_i-r n "^r.7: (2.20b) 



In the case |e| = we obtain: 

lim lOr)| = l. (2.21) 

\r\— >oo 

Relations ()2.2())1 and (|2.21j) justify the exact behavior of local solutions (|2.14|) . when 
one approaches the infinite singular point from different directions in the complex plain 
C r . They give precise meaning of limits ()2.18|) for X = oo. These relations will play an 
important role in our novel approach to different spectra of the corresponding two-points 
boundary problems related to the Regge- Wheeler equation. 



3 Two-Point Boundary Problems 

3.1 General Description and Symbols of Two-Point Boundary 
Problems for Regge- Wheeler Equation 

The two-singular-point boundary problems can be defined on the real intervals which 
contain only two singular end-points and have no other singular points inside. It is 
possible to consider regular-singular two-point boundary problems too. In this case, one 
end of the interval is a singular point and the another is a regular point supplied with 
some standard boundary condition (for example - with Dirichlet's boundary condition, see 
below). Different boundary conditions correspond to different physical problems described 
by the same differential equation. 

One can find the general definition of a stationary two-singular-point boundary prob- 
lem in the book by Slavyanov and Lay in [Sj. For our problem we need some natural 



generalization of this definition which includes not only the case of solutions, which are 
finite at both singular points. Following the physical traditions (see articles j^), we use 
as a basis of such a definition the in-out properties of the time- dependent solutions of the 
Regge- Wheeler equation, considered in the previous Section. Since we have establish the 
one-to-one correspondence between these properties and the corresponding space limits 
(12.18)1 on the real axes lR r , we can reformulate the definition using only the solutions of 
the stationary Regge- Wheeler equation ()1.2|) in the spirit of the book by Slavyanov and 
Lay. 

All two-points boundary problems can be formulated using proper constraints on the 
transition coefficients Ty±(s, s,l) in relations (|2.1fi)l . This way one obtains equations for 
finding of corresponding discrete spectra of (possibly complex) values of frequencies e. 

3.2 Two-Singular-Point Problems on the interval (l,oo) 

This interval corresponds to the outer domain of Schwarzschild black holes, where the 
eventual observer lives. Therefore, it is most important from a physical point of view 
and, correspondingly, most well studied at present. Nevertheless, even for this interval 
much further work remains to be done. 

3.2.1 Quasi-Normal Modes 

We start the illustration of the general statements, made in previous Section, with the 
most well studied at present example of quasi-normal modes (QNM) of Schwarzschild 
black holes |3], as a specific two-singular-point boundary problem on the interval (1, oo). 
Using the graphical scheme, shown in Table Q we can symbolically write down relation 
(12. 16(1 in terms of time dependent solutions of the Regge- Wheeler equation as follows: 

{1 <-} = ri+(e, s, /){-> oo} + r£+(e, s, Q{<- oo}. (3.1) 

By definition, for QNM we have no outgoing waves from the infinite point [3], i.e. 

rii(e,s,0 = =>e = e(n,s,l), n = 0,1,2,.... (3.2) 

This equation generates the spectrum of QNM. For the above values of e Eq. (J2.16)) 
reads: 

R%, sA s,i{r) = ri+( e (n, s, I), s, l)R^Jr). (3.3) 

We can simplify significantly the notation if we introduce the symbol " 1 oo" for 
this two-point boundary problem. It reflects symbolically the properties ()2.18j) of its 
solutions which enter into relation ()3.3|) . Accordingly, we denote by $iwoo(*, r\n, s, 1) = 
e * £ i-<x>(n,s,0^ lwoo (r|n, s, I) and Ri^ooi^n, s, I) = R 3^ s n s t (r) the solutions of this bound- 
ary problem, and by £i-_oo( ra , s, I), the corresponding eigenvalues. 

3.2.2 Left Mixed Modes 

Under the condition 

r£ + (e, s, Z) = =>e = ei„oo(n,s,Z), n = 0,1,2,...; (3.4) 



we obtain the spectrum of the waves, going left from infinity to horizon r = 1 without 
reflection - left mixed modes (LMM) 5 . Due to Eq. (j2.18|) we use for them the symbol 
l^oo. 

Considering the equation 

{i -} = r^-(e, s, 00} + ri:( £ , s, i){<- 00} (3.5) 

we can introduce two more types of modes of Schwarzschild black holes: 

3.2.3 Normal Modes 

The equation 

T^~(e, 8,1) = =^e = £1^00(71,3,/), n = 0,l,2,...; (3.6) 

defines the spectrum of normal modes (NM) on the interval (1, oo), which are important 
for considerations of the stability of Schwarzschild black holes. For obvious reasons, 
related once more to Eq. (|2.18j) . we use for them the symbol l^oo. 

3.2.4 Right Mixed Modes 

Finally, the equation 

r^:(£,s,0=0 =>e = e 1 „ OB (n,s,l), n = 0,1,2,...; (3.7) 

defines the spectrum of the waves, going right from the horizon r = 1 to infinity without 
reflection - right mixed modes (RMM). Due to Eq. (|2.18j) we use for them the symbol 
1 ~oo. 

3.3 Two-Singular-Point Problems on the intervals (X, Y) 

The generalization of the previous consideration for any real interval (X, Y) with singular 
ends X < Y and without other singular points in it is obvious. One can introduce the 
corresponding: 

a) Quasi-normal modes: ^x^Y(t,r\n, s,l) = e t£x ^' Y( - n ' s '^ t Rx^Y(r\n, s, I) - the solu- 
tions of the boundary problem with eigenvalues Sx-^ y{n, s, I) defined by the equation: 

r£+(£,s,Z) = =>- £ = e x ^Y(n,s,l), n = 0,1,2,.... (3.8) 

b) Left mixed modes: $x-y (t, r\n, s, I) = e l6x ^ Y( - n ' s '^ t Rx^,Y( r \ n , s, I) - the solutions 
with eigenvalues ex^y( n , s, I) defined by the equation: 

r£f(e,s,0 = ^£ = e x ~Y(n,s,l), n = 0,1,2,.... (3.9) 

c) Normal Modes: y(t, r\n, s, I) = e t£x ~ Y ( n ' s ' lS)t Rx^,Y(r\n, s, I) - the solutions with 
eigenvalues y(n, s, I) defined by the equation: 

r£-(e,s,Z) = =>e = e x ~Y{n,s,l), n = 0,1,2,... . (3.10) 

5 We use the term "mixed modes" for indication of the following obvious property of these solutions: 
Their boundary conditions at the two end points of the corresponding interval are of different type, in 
the sense of Eq. (|2.18|l . 



d) Right mixed modes: $x~y (t, r|n, s, I) = e l£x ~ Y(jL ' s, ^ t Rx~Y{f\ n ^ s ? ~~ t^e solutions 
with eigenvalues Ex~y( n , s, I), defined by equation: 

rJl(e,s,Z) = =>-e = ex~y(n,a,0, n = 0,1,2,.... (3.11) 

To the best of our knowledge, all spectra of Schwarzschild black hole two-point bound- 
ary problems on the intervals (X, Y) ^ (1, oo) at present are unknown both analytically 
and numerically. Much better is the situation in the case of the interval (l,oo): the 
spectrum of QNM is known in detail, both numerically and analytically 3J, although no 
exact explicit expressions have been obtained up to now. In addition, in this case we have 
some qualitative studies of the spectrum of normal modes [§]. 

3.4 Regular- Singular- Two-Point Boundary Problems 
3.4.1 Formulation of the Problem and Notations 

These two-point boundary problems for the Regge- Wheeler equation (jl.ljl are related to 
models of heavy compact objects which are essentially different from a black hole model. 
The problem of finding the corresponding spectra in the simple mathematical formulation, 
considered below, seems to be new and even not posed in the existing physical literature. 

In the present section we consider a one-parameter family of real intervals (r*, oo) C 
(1, oo). The left end of the interval is an arbitrary regular point r* G (1, oo) of Eq. (jl.2|) . 
We pose Dirichlet's boundary condition at this end: 

r.) = 0. (3.12) 

Physically, this means that the waves are experiencing a total reflection at the spherical 
surface with the area radius r = r* > 1, instead of going freely through it without any 
reflection, as in the case of black holes. We have to mention that Dirichlet's boundary 
condition at the cosmological horizon has been used recently for QNM of AdS black holes 
in the article jS]. 

One can consider this surface as a boundary of some massive spherically symmetric 
body. Owing to the Birkhoff theorem, in general relativity the gravitational field of any 
massive body with Keplerian mass M outside the radius r* coincides with the correspond- 
ing field of a black hole. Hence, the spreading of the waves in the outer domain of the 
massive body is governed by the Regge- Wheeler equation fjl.ip . too, but the boundary 
conditions must be different. The corresponding effective potential V Sj i(r) in the Regge- 
Wheeler equation (jl.ljl is illustrated in Fig. |2J 

There exists a large amount of articles on QNM of relativistic stars [31 Ej, both sta- 
tionary and rotating. At present, this field is well studied in detail. It is known that, in 
general, the interior of the stars may play a complicated role in the polar oscillations of 
fields of different spin outside the star. Some details are still an open problem, because 
the relativistic equation of state of matter is not precisely known for the real case of very 
high pressures in the compact relativistic stars, for example, in neutron stars. 

For axial perturbations of metric the situation is different: axial perturbations do not 
induce a motion of the star matter and experience only a potential scattering, as in the 
case of Schwarzschild black holes [3]. This is a physical reason for our simple model of 
axial perturbations of gravitational field of compact objects. 
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Figure 2: The potential V s j(r) for s = 2, I — 2. The small circle marks the place of 
its maximum (here and at all figures hereafter). The shadowed domain marks the area, 
which is forbidden for spreading of waves. Its right edge is at the point r*. The point 
r = 1 presents the event horizon. 



In this model, we avoid all complicated problems related to the body's interior, replac- 
ing them with the simple boundary condition (J3.12j) . It can be thought as idealization, 
which describes a proper limiting case of problems related to real bodies. It is obvious 
that this idealization is opposite to the Schwarzschild black hole model. If there ex- 
ist a conserved current in the problem at hand, one can expect that the corresponding 
eigenvalues for real static spherically symmetric objects will be placed between the cor- 
responding eigenvalues of black holes and the eigenvalues of the simple model of massive 
compact objects, we study here for the first time. Unfortunately, such simple guess turns 
to be incorrect, because the current in the RW problem with complex eigenvalues e is not 
conserved and has quite nontrivial properties, see the Appendix. 

In general, one can try to consider two different boundary problems: 

a) Quasi-normal modes of massive compact objects with the spectrum e^^^n, s, I): 
$n~oc{t, r\n, s, I) = e* 6 * --fr^i^ *,oo(r|n, s, I), n = 0, 1, 2, . . . . 

b) Normal modes of massive compact objects with the spectrum e n »^oo(n, s, I): 
$ n *oo(t,r\n,s,l) = e u ^- oo ^ n ' a ' t > t R umn0O {r\n > s,l), n = 0, 1,2, . . . . 

Here we use the symbols r* »~oo and r* •^oo for QNM and NM in the exterior domain 
of compact massive objects. In these symbols the sign "r**" denotes the spherically 
symmetric compact object with the boundary of the area radius r*. The tail denotes the 
corresponding perturbation mode in accord with the properties described by Eq. (j2.18|) 6 . 

Using the properties of the complex current in RW problem one can prove that NM 
do not exist under Dirichlet's boundary condition ()3.12|) . see the Appendix. This is well 
understandable from physical point of view, since in the problem at hand we have no 
absorbtion of the waves, but a total reflection, creation and emission to space infinity. 

6 From a pure mathematical point of view it is possible to consider analogous modes on intervals (r + , X) 
with X = l,r* G (0, 1) and X = l,r* S (l,oo), orl = 0,r, £ (0, 1) and even X = 0,r* G (-oo, 0), using 
symbols like r* »^X and r* »~X for the corresponding QNM and NM. Since the physical meaning of such 
modes is not clear, we will ignore them in the present article. 



3.4.2 The Equation for the Discrete Spectrum of QNM 



Let us consider solutions of the Regge- Wheeler equation <3> £ r * r) on the interval (r*, oo), 
which obey the boundary condition (j3.12j) . Making use of solutions (j2.13b|) of the station- 
ary equation (jl.2|) . we can represent the solutions $ £ r * ; (t, r) in the following determinant- 
form: 

= ^W$SJ(*' r )-^i"!i(r.)*S,T(*' r )- (3-13) 

Their expansion on the basis of solutions (|2.14|1 . represented in the symbolic form 
reads: 

r) = r* + (e, a, 00} + r s, /){- oo}, (3.14) 

where r o ^* ± (e, s, Z) are proper coefficients. It is obvious that the spectrum of QNM is 
defined by the following equation: 

r*_(e,s,Z) = e ft ^oo(n, 3,1), n = 0,1,2,.... (3.15) 

Now one can use relations (|2.15j) and (|2.1fjj) to obtain: 

T^(e, s, I) = T^+(s, s, l)R%(n) - T^(e, s, l)R { l ] +(n) (3.16) 

and to represent Eq. ()3.15|) in the form: 

r^M,z) _ iQn) _ 

c2l£Mn . 1) HeunC(-2ie,+2ie,2s-2e 2 ,2e 2 + s 2 -l(l + l),l-n) 

HeunC(-2ie-2ie,2s-2e 2 ,2e 2 + s 2 -l(l + l),l-r*)' 1 ' ' 

Here we have used the relations (|2.13b|) . too. Relation ()3.17|) present the final form of the 
equation for the spectrum of QNM in our simple model of massive compact bodies. 



R%ir) 



3.4.3 Some General Properties of Trajectories of QNM Eigenvalues 
£r* «-oo(^, s, Z) in the Complex Plane C^ ^^ 

For the massive compact bodies with different values of the area radius of their surface 
the eigenvalue e n ^^(n, s, Z), as a function of G (l,oo), runs on some trajectory in 
the complex plane C £rfc#iioo . Here we prove three simple properties of these trajectories. 
These properties are confirmed in the next section by numerical analysis of trajectories 
and illustrated in Figs. ® and (jH). 

Proposition 1: If the limit r* — > 1 + of the trajectory of e^^^n, s,l) in the complex 
plane C r exist, then we have: 

lim e n ^oo{n, s, I) = 0. (3.18) 

r*— »l+0 

Proof: The eigenvalue e n ^oo( n , s, Z) is obtained from the corresponding Eq. f)3 . 1 7j) . 
Hence, its dependence on r* is defined by the rhs of this equation. 



The lhs of Eq. ()3.17j) must be nonzero and finite, because the solutions ()2.13j) and 
(12. 14(1 are two complete sets of independent local solutions of the Regge- Wheeler equation. 

Suppose that the limit (j3.18|) exists. According to the relations (|2.9aj) and (j2.13b|) . 
the limit r* — > 1 + of the rhs of Eq. (|3.17j) : 

does not exist if lim^^i+o (^(sr,^oo( n , s > 0)) 0- If h m »wi+o ffi{ £ n «-oo( n > s , 0)) = 
and lim rt ^ 1+0 ($s(£ n »~oo( n , s, I))) > 0, the limit ()3.19|) is infinite. In the case 
lim^-si-H) (S*(£r» ^oa(n, s, I))) < the limit ()3.19j) is zero. Hence, if a nonzero finite limit 
(l3~TT3j) exists, Eq. (jXTSjl must be fulfilled.^ 

Proposition 2: In the limit r* — ► +oo on the trajectory of s^^^n, s, I) in the 
complex plane C eii#JOO we have: 

lim e n ^oo(n, s,l) = 0. (3.20) 

r*— »+oo 

Proof: Using the expansion on R^f(r*) of type (|2.16|) for ^^(r,) in Eq. (j3.17j) and 
the properties ()2.20|) . one easily obtains the limit r* — > +oo of both sides of relation (j3.17j) 
in the form: 

lim f^<*>-( £r * S ' 0' S 5 \ _ ^ m /I 1 oo+( £ r* •^oo( ri ) s ; 0) S ' A (3 21) 

r,^+oo\T^ Z(e n ^ OD (n,s,l),s,l)J *-h-og Vr^(e r ^ 00 (n, s, I), s, I) J 

This means that in the limit r* — > +oo the solutions R^^r) become linearly depen- 
dent. Formulae ()2.13|) show that this is possible only when equation ()3.20j) takes place, 
o 

Proposition 3: The trajectory of eigenvalue e n ^oo(^, s, I) in the complex plane 
C £rat#JOO has no cusps, i.e. the real part 3? (e n ^oo(n, s, I)) and the imaginary part 
s,l)) can not have simultaneously extremum for finite values of r* G (l,oo). 
Proof: Let us introduce the short notation g(e) = T^(e, s,l)/T^ >D Z(s, s,l) and 

h(s,r*) = R^*i(r*) I R^lsii 1 "*) • Then Eq. (I3.17|) . written down in the short form 
g(e) — h(e, r*) = 0, defines implicitly the function sir*) in the general case when d e ( K g{e) — 
/i(£,r*)) 7^ 0. For a derivative of the function e(r„) one obtains d £ (g(e) — /i(e, r*)) ^ = 

dh(e,n) = -W (R^(n),R%(n)) / (<^(r,)) 2 , where W (it^r*), R%(n)) de- 
notes the corresponding Wronskian. Since ^^s~X r *)) 7^ oo for finite r* G (1, oo), we 
can have ^(f*) = if, and only if W f-R^jO*)) R^^i( r *)) — 0. This contradicts to linear 
independence of the local solutions R^lsii 7 *) and i?^ s ^r*) for finite r* G (1, oo). o 

4 Numerical Calculation of Quasi-Normal Modes of 
Static Spherically Symmetric Objects 

4.1 Transition Coefficients From Limit Procedures 



As we have seen in the previous sections, for determination of QNM in simple models of 
compact objects we need to use the transition coefficients Ty±(£, s,l). Since these are 



not known explicitly, one has to apply more sophisticated methods for solving two-point 
boundary problems, instead of employing directly Eqs. (|3.2|) and (j3.17l) . 

In this section we present a new simple technique for calculation of quasi-normal modes 
of Schwarzschild black holes and massive spherically symmetric bodies, using directly the 
exact solutions of the Regge- Wheeler equation, found for the first time in the present 
article. We use our own computer code, written in Maple 10 package. There a pioneering 
implementation in computer algebra systems of the five Heun functions for complex values 
of argument and of corresponding parameters already exist. We use the numerically 
well studied case of Schwarzschild black holes to compare our results with the earlier 
ones, which were obtained by making use of two other numerical methods of Leaver and 
Andersson 

The basic idea of our new method is to use proper limit procedures for calculation of 
transition coefficients. 

Let us consider some solution R EjS} i(r) of the stationary Regge- Wheeler equation and 

let 

R^r) = C 1+ R%(r) + C^R^r) = CU<i>) + C^R^J(r) (4.1) 
present its expansions on the corresponding basic solutions fl2.13j) and (|2.14|) . Then 

Oo± = C 1+ T^ ± (e, s, I) + C^T^e, s, /), (4.2) 

and we obtain the modulus of these constants using relations ()2.2())) in calculation of the 
following limits: 

ICUI = lim \e + ^ r+la ^R £ ^(r)\ : for arg(r) + axge G (0, n), (4.3a) 

\r\— >+oo 

lOo^l = lim \e- ie{r - lnr) R £ , Stl (r)\: for arg(r) + arge G (-tt, 0). (4.3b) 

\r\— >+oo 

In addition, it becomes clear from asymptotics ()2.19|) that in the limit (j4.3b|) we have 
the steepest descent, if and only if arg(r) + arge = — | / mod 2ir. This direction is the 
optimal one for numerical calculations of the above limits, as we shall see below. 

This way we can obtain the needed information about transition coefficients 
T ^±(e , s , I) , using any solution with known constants C\±. Alternatively, we can im- 
pose the corresponding boundary conditions directly on the modulus |Coo±|, calculated 
via the limits (|4.3j) . 

4.2 QNM of Schwarzschild Black Holes 

Already in a first calculation the QNM frequencies, made by Chandraseckhar and De- 
tweiler [Z| the following difficulties were recognized: It is hard to nullify numerically the 
coefficient in front of the solution R^si ( r ) ( see ^q- ()3.2j) ). which is small on the left of 
the infinite point on the real axes. This happens because it is hard to distinguished the 
small solution R^~[(r) in the presence of unavoidable numerical errors in the solution 

R^Tst( r )> w hich is dominant in this domain. Therefore, much more sophisticated meth- 
ods, like Leaver's continuous fraction method and Andersson's phase amplitude method 
for accurate calculations of QNM spectrum, were invented successfully, see references in 
[7j. This way were obtained the known most accurate values for QNM frequencies. 



For example, in Andersson's article, among many others, one can find the value of the 
basic quasi-normal mode of gravitational waves: £i^oo(0, 2, 2) Anders son = 0.747343368 + 
% 0.17792463. 

Using the Leaver's method on modern computers, it is easy to find several more figures: 
£i^oo(0,2, 2) \ Leaver = 0.747343368836 + i 0.177924631377 for the same mode 7 . 

Since today we have at our disposal a possibility to calculate the canonical solution 
(J2.3J) of Heun's equation in the complex domain of variables using Maple 10, we are 
able to turn back to the idea of solving numerically Eq. (|3.2Jl ). in a properly modified 
way, which overcomes the difficulty of Chandraseckhar-Detweiler method. Our solution is 
based on the limit procedure (J4.3bj) which shows that in the direction arg{r) + arg£ = — § 
the solutions interchange their roles. The former small solution becomes a dominant 
one. Taking into account expressions (|2.13|) . one easily obtains the following form of the 
equation for for QNM spectrum of Schwarzschild black holes: 



lim 

|r|^+oo 



s+1 HeunC (-2ie, 2ie, 2s, -2e 2 , 2e 2 + s 2 -1(1 + 1), 1- |r|e<^ rg{e) ) 



= 0. (4.4) 



In a numerical treatment of the problem we replace the last limiting procedure with 
calculation at some finite point ^> 1, i.e. we use the equation: 



) s+l HeunC (-2ie, 2%e, 2s, -2s 2 , 2e 2 + s 2 - 1(1 + 1), 1 - r 00 e<^ +arg{£) ) 



0. (4.5) 



We apply the complex plot abilities of Maple 10 to find a relatively small vicinity of the 
complex roots of this equation. Then we justify their values to the desired precision using 
a direct check of the place of the corresponding root in the complex plane C e and/or by 
Muller's method for finding of complex roots 10 . 

The results are shown in the Tables El El and El In the Table El the results of our 
numerical calculations for the first five eigenvalues are presented. We see an excellent 
agreement with the most accurate values, published up to now by Andersson jjj. 



n = 


0.747343368 + i 0.177924631 


n — 1 


0.693421994 + i 0.547829750 


n = 2 


0.602106909 + i 0.956553966 


n = 3 


0.503009924 + i 1.410296405 


n = 4 


0.415029159 + i 1.893689781 



Table 2: The first five eigenvalue for 64 digit calculations and for = 20. 



The Table El shows that using high precision calculations with Eq. ()4.5)1 one can 
produce the value of £1^00(0, 2, 2) up to more figures - here up to 12 figures. This 

7 The author is deeply indebted to prof. Kostas Kokkotas for permission to use his computer code for 
calculation of this value by Leaver's method. 



consideration inspires confidence in our new method for calculation of QNM. As seen 
from Table E3 with our Maple 10 computer code we obtain the right 9 figures already 
in 64 digits calculations. This gives us more confidence in the results for next several 
eigenvalues, obtained using the same method and shown in Table [21 



32 


0.74734336c 


3602 


+ i 0.177924631717 


64 


0.74734336c 


3782 


+ i 0.177924631969 


128 


0.74734336c 


3864 


+ i 0.177924631732 


256 


0.74734336c 


3836 


+ i 0.177924631360 



Table 3: The basic eigenvalue for higher number of digits (given in the first column) 
calculations (no rounding) and for = 20. 

In the TableEJwe show the deviations of the values of the basic eigenvalue £i^oo(0, 2, 2), 
obtained by our method for different r^, with respect to the corresponding value for 
= 100 and for 32 digits calculations. 



10 


-.17e-8 - i .211e-7 


100 


+ i 


20 


-.90e-8 - i .106e-7 


200 


.252e-7 - i .141e-7 


30 


-.96e-8 - i .105e-7 


300 


-.156e-7 + i .204e-7 


40 


-.65e-8 - i .82e-8 


500 


.2007e-6 - i .3180e-6 


50 


-.103e-7- i .160e-7 


1000 


-.785e-7 - i .3198e-6 


80 


-.94e-8 - i .50e-8 


1500 


.2680e-6 - i .21134e-5 



Table 4: The numerical results for 32 digit calculations of £i^oo(0, 2, 2) for different values 
of ^ (given in the first and in the third column). In this Table we present only the 
deviations from the value e^^O, 2, 2)| 100 = 0.747343377577 + i 0.177924642293, 
obtained for = 100. 

As we see, the calculation of the eigenvalue £1^00(0,2,2), based on equation (|4.5|) . is 
not very sensitive to the choice of and for > 20 it gives a stable right result up 
to seven figures for 32 digits calculations. A simple look on the Table 0] shows that the 
deviations have no systematic nature and may demonstrate instabilities (after the first 
seven figures) of Maple 10 calculations with 32 digits precision. 



4.3 QNM of Massive Compact Body 



Applying the limit procedure (j4.3b|) to our simple model of massive compact body, de- 
scribed in Section 3.4, we obtain the following equation for the corresponding QNM: 



lim 1 < 

Irl— >+oo 



2i*e ln(|r|e- i (' r / 2 + ar s( e ))-l) 



X 



HeunC (-2ie, +2ie, 2s, -2e 2 , 2e 2 + s 2 - 1(1+1), 1- 


r 


e -i(n/2Urg{e)J\ 


HeunC {-2ie, -2ie, 2s, -2e 2 , 2e 2 + s 2 -l(l+l), 1 - 


r 


e -i(7r/2+arg( £ ))) 



HeunC (-2ie, +2ie, 2s, -2e 2 , 2e 2 + s 2 - 1(1+1),!-^ 



(4.6) 



e 2i*eln(r*-l) x 



HeunC (-2ie, -2ie, 2s, -2s 2 , 2e 2 + s 2 -1(1 + 1), l-< 



0. 



The numerical treatment of this problem is the same as in the previous Subsection: 
we once more replace the limiting procedure with calculation at some finite point ^> 1, 
i.e. we use the equation: 



_2i*e ln(r 00 e-*<> r / 2 + ar st e ))-l) 



HeunC (-2ie, +2ie, 2s, -2e 2 , 2e 2 + s 2 - 1(1 + 1), l-r^e-^^ 5 ^) 
HeunC (-2ie, -2ie, 2s, -2e 2 , 2s 2 + s 2 -l(l + l),l-r O0 e- i ^/ 2 ^^ £ ») 



(4.7) 



2i*e ln(r,— 1) x 



HeunC (-2ie, +2ie, 2s, -2e 2 , 2e 2 + s 2 - 1(1 + 1), 1-r*) 



HeunC (-2ie, -2ie, 2s, -2e 2 , 2e 2 + s 2 -1(1 + 1), 1-r*) 



and apply our own computer codes, written in Maple 10, for finding its roots. The results 
are shown in Fig. Eland Fig. HJ 



5 Conclusion 

In the present article we have solved some mathematical problems related to exact so- 
lution of the Regge- Wheeler equation. Since the obtained results and possible further 
developments were described in detail in the basic text, in the concluding remarks we will 
comment on their eventual physical significance and some open problems. 

In the present article we do not establish a definite relation between the spectrum of the 
problem with Dirichlet's boundary condition at the surface of the spherically symmetric 
body and the corresponding spectrum of more realistic models of relativistic stars 0. 
The simple comparison of the published numerical results for relativistic stars with our 
results in Figs. El lU shows that the complex eigenvalues of star's axial modes are placed 
in a complicated way in the domains, separated by the trajectories of QNM eigenvalues 
£r**<oo(n, s, I), which belong to the simple model of massive compact objects, studied here 
for the first time 8 . One may hope that further investigation of this problem will help our 
understanding of QNM of compact matter objects. 

A more realistic model for real spherically symmetric bodies may be obtained consid- 
ering as a boundary condition a partial reflection and partial penetration of the waves 

8 The author is grateful to unknown referee for its suggestion to consider such comparison. 
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Figure 3: The complex trajectory of the basic eigenvalue .^(O, 2, 2) of QNM of 
compact massive body as a function of the distance d = r* — 1 > of the body's surface 
to the horizon (in units 2M). For comparison purposes, the first two eigenvalues of QNM 
for Schwarzschild black holes with the same mass M are shown, using the circles with 
crosses. The origin of the complex plane C £ is marked by a small cross. It presents 
simultaneously the event horizon (where the trajectory begins) and the physical infinity 
(where the trajectory ends). The small circle on the trajectory denotes the position of 
the maximum of the effective potential V(r). 
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Figure 4: The trajectories in the complex plane C e of the first three eigenvalues: 
£V*~oo(0, 2, 2),e r , >v , 00 (l, 2, 2) and £,.„, .^(2, 2, 2) of QNM of massive compact body. For 
comparison purposes, the first six eigenvalues of QNM for Schwarzschild black holes with 
the same mass M are shown, using the circles with crosses. The small circles on the tra- 
jectories denote the corresponding positions of the maximum of effective potential V(r). 



through the body's surface. In the present article we have considered only the two extreme 
cases for reflection coefficient: R — - for the black holes and R — 1 - for our simple 
model of bodies. It is interesting to know will it be possible to connect continuously 
these two cases, introducing a coefficient of reflection R G [0, 1], or the difference in the 
boundary conditions is an insuperable obstacle for this. 

In the near future our results may be useful for the planned observations of grav- 
itational waves emitted by or spreading around existing real compact objects. These 
results may help to clarify the physical nature of the observed very heavy and very com- 
pact dark objects in the universe. Our consideration gives a unique possibility for a 
direct experimental test of the existence of space-time holes. The study of spectra of the 
corresponding waves, propagating around compact dark objects, may give indisputable 
evidences, whether there is space-time hole inside such an invisible object. 

Similar results can be obtained for electromagnetic, spinor and scalar waves in the 
external Schwarzschild gravitational field of compact objects. A further study of these 
phenomena is an important physical problem and may offer practical physical criteria for 
experimental distinction of the two completely different hypothetical models: 

1) The model of space-time black holes; and 

2) Different new models of very compact dark objects, made of real, or of some hypo- 
thetical matter jTT] . 

The present-day astrophysical observations still are not able to make difference be- 
tween these two cases. The only real observational fact is that we see very compact and 
very massive objects which show up only due to their external strong gravitational field. 
An actual theoretical problem is to find a convincing model for description of these al- 
ready observed compact dark objects and the criteria for experimental verification of a 
model like this. We hope that the present article is a step in this direction. 
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6 Appendix: General Properties of the Complex 
Current and Some Important Consequences 

6.1 General Properties of the Complex Curl 

At first glance the Regge- Wheeler eigenvalue problem seems to be quite similar to the 
Schrodinger eigenvalue problem in quantum mechanics. There one considers usually real 
eigenvalues of the energy operator. The last assumption yields existence of a conserved 
real current. This current describes an essential physical information about the solutions 
of the problem. We were not able to find in the literature a proper investigation of the 
properties of the current in RW problem with complex eigenvalues e. Therefore in this 
Appendix we define the corresponding complex current, study its basic properties and 
derive some consequences, which are essential for our article. 

In the case of e G C, er > 0, Ej ^ 0, and 7gR the stationary RW equation (jl.2j) for 
R = X + iY = \R\e %ip is equivalent to the real system of 4th order ordinary differential 
equations: 

X" — (V — OA + V Y, Y" = (V - £)Y - rjX, (6.1) 

where the prime denotes differentiation with respect to RW tortoise variable x and e 2 = 
e 2 r — e 2 + 2ierEy = £ +irj. In the case i] ^ this system does not split into two independent 
differential equations of 2nd order and in addition we obtain: 

r _ v , (X') 2 + (Y>) 2 (XX' + YY')' _ (YX' - XY')' 

* X 2 + Y 2 X 2 + Y 2 ' V X 2 + Y 2 ' 1 ' ' 

Despite of the above fact, for r\ ^ the system (|6.1|) is equivalent to 2nd order system 
of real nonlinear differential equations (see, for example, the article by S. Chandrasekhar 
and S. L. Detweiler in [7j): 



u 



R 



7] - 2u R u u u'j = u\ - u\ + V - f , (6.3) 



with boundary conditions at infinity: ur(oo) = Er > 0, 1*7(00) = Ej > - for QNM. 
This phenomenon turns to be possible since two of the arbitrary real constants in the 
general solution of the system (j6.1|) are actually not essential, entering inessential com- 
plex constant multiplier of form i?* = \R*\e llp * in the solution R of the RW eq. ([1.2)1 . 
This constant always can be chosen to be equal to 1. To obtain Eqs. ()6.3|) we use the 
substitution: 

R = e ~iju{x)dx _ e -i J u R {x)dx e J Ul (x)dx ^ (QA) 



where 

.R' .R'R jR + tji 
u = u R + i Ul = i—=i—= . (6.5) 

This way we introduce a complex current j = j R + ijj = u R \R\ 2 + iui\R\ 2 in the RW 
problem. Using the Eqs. (jfi.lj) . it is easy to obtain the following relations: 

j R = YX'-XY', f R = n(X 2 + Y 2 ), j' R = 2 V (XX' + YY'), 

ft = 277 ((X') 2 + (F') 2 ) + 2 V (V - 0(X 2 + F 2 ), (6.6) 

and the 4th order differential equation for the real part j R of the complex current 

Jr V = <V - t)j" R + V% + A V 2 JR . (6.7) 

One sees that for n ^ 0: 

i) The complex current j is not conserved, i.e. j R ^ const and ji ^ const. 

ii) Since jj = 2~j R = XX' + YY', in the imaginary part jj of of the current j there is 
no independent information. 

hi) In addition \R\ 2 = j' R /n, tp' = —nj R /j' R . Hence, the real component j R of the 
complex current j contains all essential information for the RW problem and the Eq. 
(|6.7|) can replace RW one (jl.2|) in all considerations. 

It is easy to obtain the first integral \ {{j' R ) 2 )" — i{j R ) 2 — {V — £){j' R ) 2 — n 2 j R = const 
of the Eq. (|6.7|) . The constant value of this integral turns to be one of the two inessential 
constants in the general solution of the Eq. (J6.7|) . Without loss of generality we can 
assume that this constant equals zero. Then the real part j R obeys 3th order differential 
equation: 

\mr = l^f + (V-Wn) 2 + V 2 f R (6.8) 

Using the substitution j R = exp (n J u R (x)dx) one can exclude the second inessential 
constant in the general solution of the Eq. ()6.7|) . lowering once more its order. Thus we 
obtain the single 2nd order differential equation: 



1 fv ~ u' R \ 1 { v - u' R 



2 



, » , +«£ + V-£, (6.9) 

2 V U R J 4 V u R J 

which can be derived directly from the system ()6.3J) . too. The last observation supports 
our choice of the zero value for the first integral of Eq. ()6.7|) . 

At the end let us introduce functions R^ — | [R + jlR') , R^ = \ [R — hM) putting 

R = R^ + R_>, R' = ik(R^-R^). (6.10) 



Here k = ye 2 — V and we suppose to choose the branch of the square root for which we 
have k = e when V = 0. Now we obtain 

j = fc(|i^| 2 -|i^| 2 ), (6.11) 

precisely as in the real case ej = 0. 

The complex functions R^ and R^. obey the 2nd order system 

£)'-*(£)■ M1!)*m')' (612) 

which shows that at the points where V — these functions describe waves which run 
to the right and to the left, correspondingly. If, in addition, V — 0, then = e R \R^\ 2 , 
= -e R \R^\ 2 . 



6.2 Special Properties of the Current in the Case of Compact 
Body with Dirichlet's Boundary Condition at the Surface 



In the simple model of compact body with Dirichlet's boundary condition (|H.12j) one has 
R(x*) = 0, R'x = -R'( x *) 7^ 0- As a result in vicinity of the point x* one obtains the 
following series expansions with respect to the difference Ax = x — x*: 



R(x) = R\ {ax + l(V*~ s 2 )Ax 3 + ^Ax 4 + 6 (Ax)^) 



(6.13) 



u R {x) = ~Ax + 3 (Ax), (6.14a) 

Ul {x) = J- + i(K - i)Ax + ]v; Ax 2 + 3 (Ax); (6.14b) 
Ax 6 4 



jR (x) = V \K\ 2 ( ^Ax 3 + - Ax 5 + ^K'Ax 6 + 7 (Ax) ) , (6.15a) 



j 7 (x) = K| 2 (^Ax + -(K - 0Ax d + -K'Ax 4 + 5 (Ax) ) ; (6.101,) 
and the two equivalent sets of integral representations: 

X 



u R (x) = r]\R(x)\- 2 j \R(x)\ 2 dx, (6.16a) 

X* 
X 

m(x) = \R{x)\- 1 J {\u R (x)\ 2 + {V(x) - t)\R(x)\) dx + (6.16b) 

and 

a; 

3R{x)=r] J \R(x)\ 2 dx, (6.17a) 

X X 

3l {x) = J (\R'(x)\ 2 + V{x)\R{x)\ 2 )dx-i J \R{x)\ 2 dx. (6.17b) 
With the help of the Eq. (I6.17al) the Eq. (I6.17bl) can be rewritten in the form 

X 

&r{x) + VJi(x) = V ( {\R\x)\ 2 + V(x)\R(x)\ 2 ) dx. (6.18) 



Now we obtain the following additional properties of the solutions of RW equation for 
compact body with Dirichlet's boundary condition at its surface: 



1. NM do not exist for rj 7^ in the Dirichlet's boundary problem ()3.12j) . Indeed, 

00 

since for NM R(oo) = 0, one must have for them = jr(oo) — rj J \R(x)\ 2 dx. Hence, 

Xt 

R(x) = for x G [a?*, 00) and we will have a trivial zero-solution. 

2. For QNM we have r\ > 0, because of our assumption Er > 0. Then we obtain the 
following properties of QNM in this problem: 

i) The solutions R(x) and the current Jr{x) cannot have other zeros xq, different from 

x 

x*. Indeed, if R(x )—0, then = jn(x ) = r\ J \R(x)\ 2 dx > leads to x = x* for any 

nontrivial solution R(x) 7^ 0. 

ii) If Si > er, then £ < and Eq. (j6.17b|) shows that ji(x) > and j" > for any 

x e (x*, 00). 

hi) Oscillations of ji(x) (and j") around the zero value: 

From Eq. (j6.17b|) we see, that ji(x) (and j R (x) = 2r]j I (x)) can have a second 
zero, besides x„, if and only if Ej < Er, i.e. when £ > 0. Then, since lim ]r(x) = 

x— >oo 

lim (ur(x)\R(x)\ 2 ) = +00, the current Jr{x) must have at list one more (third) zero in 

X—* OO 

the interval [x*, 00). This way we see that in general case for £ > the imaginary part 
ji(x) of the complex current has an odd number (> 1) of zeros and its derivative jj{x) 
has an even number (> 0) of zeros in the interval [ 

iv) From Eq. (j6.18|) for x > x* we obtain the inequality: 

5(^-^7) ( 6 - 19 ) 
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